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ABSTRACT 

. 

The axial mode oscillations are examined for relativistic rotating stars with 
uniform angular velocity. Using the slow rotation formalism and the Cowling 
approximation, we have derived the equations governing the r-mode oscillations 
up to the second order with respect to the rotation. In the lowest order, the 
allowed range of the frequencies is determined, but corresponding spatial function 
is arbitrary. The spatial function can be decomposed in non-barotropic region 
by a set of functions associated with the differential equation of the second- 



> 

in 



order corrections. The equation however becomes singular in barotropic region, 
^ . and a single function can be selected to describe the spatial perturbation of the 

lowest order. The frame dragging effect among the relativistic effects may be 
significant, as it results in rather broad spectrum of the r-mode frequency unlike 
in the Newtonian first-order calculation. 
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c/3 ■ 1. Introduction 

In recent years, the r-mode oscillations in rotating stars have been found to be signifi- 
cant implications. The axial oscillations are unstable by the gravitational radiation reaction 
( |Anclcrsson 1998] , [Friedman and Morsink 1998 ). The mechanism of the instability can be un- 



derstood by the generic argument for the gravitational radiation driven instability, so-called 
CFS-instability, which is originally examined for the polar perturbations, ( |Chandrasekhar 



1970 . Friedman and Schutz 1978] , Friedman 1978 ). All rotating stars become unstable in 



the absence of internal fluid dissipation, irrespective of the parity modes. Viscosity however 
damps out the oscillations, and stabilizes them in general. The polar f-mode instability is 
believed to act only in rapidly rotating neutron stars flLindblom 1995|) . The axial instability 
is however found to set in even in much more slowly rotating cases, and to play an important 
role on the evolution of hot newly-born neutron stars ( |Lindblom et al. 1998| , [Andersson et 



al., 1998[) . The instability carries away most of angular momentum and rotational energy 
of the stars by the gravitational radiation. The gravitational wave emitted during the spin- 
down process is expected to be one of the promising sources for the laser interferometer 



gravitational wave detectors ( |Owen et al. 1998]) . 
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Most of the estimates for the r-mode instability are based on the Newtonian calculations. 
That is, the oscillation frequencies are determined by inviscid hydrodynamics under the 
Newtonian gravity, and the gravitational radiation reaction is incorporated by evaluating the 
(current) multipole moments. Relativity has great influence on stellar structures, redshift 
in oscillation frequencies, frame dragging, radiations and so on. It is important to explore 
the relativistic effects, which may not change the general features of the oscillations, but 
shift the critical angular velocity. The relativistic calculation is not so easy task, even using 
the linear perturbation method concerning the oscillation amplitude. Both relativity and 
rotation complicate the problems considerably. 

As the first step toward a clear understanding of relativistic corrections, we have ex- 
amined the r-mode oscillations with two approximations, the slow rotation approximation 
([Hartle 1967[) and the Cowling approximation (|Cowling 1941[). The rotation is assumed to 



be slow, and treated as small perturbation from the non-rotating case. We also assume that 
the perturbation of the gravity can be neglected in the oscillations. The accuracy of the 
Cowling approximation are tested in the rotating relativistic system ( [Finn 1988] , |Lindblom| 
and Splinter 1990 , Yoshida and Kojima 1997|) . The calculations give the same qualitative 



results and reasonable accuracy of the oscillation frequencies, as in the Newtonian stellar 
pulsation theory (e.g., |Cox 1980| ). In Section 2, we summarize the basic equations for the 
inviscid fluid in slowly rotating relativistic stars and the perturbation scheme to solve them. 
The r-mode solution is given at the lowest order with respect to the rotational parameter 
in Section 3. The mode can not be specified by a single frequency unlike in the Newto- 
nian r-mode oscillation, although the frequency is bounded to a certain range. The spatial 
function of the mode is arbitrary at this order. In order to determine the radial structure, 
we need the rotational corrections up to the third order in the background. In Section 4, 
we include the rotational effects, and derive the equation governing the r-mode oscillations 
in non-barotropic region. The resultant equation becomes singular in the barotropic case, 
since the second-order differential term vanishes. In Section 5, we separately consider for 
the barotropic stars the mode function and the correction to the leading order. In Section 
6, we discuss the implication of our results. We use the geometrical units of c = G = 1. 



2. Pulsation equations of axial mode 

We consider a slowly rotating star with a uniform angular velocity Q ~ 0(e). The 
rotational effects can be treated as the corrections to the non-rotating spherical star. We 
will take account of the corrections up to the third order in e. The metric tensor for describing 
the stationary axisymmetric star is given by flHartle 1967|JChandrasekhar and Miller 1974 , 
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Quintana 1976|) 



ds 2 



-e u [1 + 2(h + h 2 P 2 )} dt 2 + e A 
+ r 2 (l + 2k 2 P 2 ) \d9 2 + sm 2 e 



2e A , 

1 H {mo + m 2 P 2 ) 
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- [u + Wi-Wa 



dr 2 

1 dP 3 " 
sin 6 d6 



dt 



where Pi(cos8) (I = 2, 3) denotes the Legendre's polynomial of degree /. The metric functions 
except g tt f) should be expanded by an even power of e, while g t< j, by an odd power of e due to 
the rotational symmetry, i.e., ui ~ 0(e), (ho, h 2 , m , m 2 , k 2 ) ~ 0(e 2 ) and (W\, W3) ~ 0(e 3 ). 
These quantities are functions of the radial coordinate r. 

The equilibrium state is assumed to be described with perfect-fluid stress-energy tensor. 
The 4-velocity of the fluid element inside the star has the components, which are correct to 
0(e 3 ), 

U+ = nU\ U l = (-g tt - 2Qg t4> - Q 2 g H )- 1/2 . (2) 

The pressure and density distributions are subject to the centrifugal deformation, which is 
the effect of 0(e 2 ). These distributions are expressed as 

V = Po(r) + {p2o(r) + P22(r)P2(cos9)}, (3) 

P = Pair) + {p2o(r) + p 22 (r)P 2 (cos#)}, (4) 

where po an d Po are the values for the non-rotating star, and the quantities in the braces 
are the rotational corrections. The non-rotating spherical configuration and the rotational 
corrections in eqs. ([[])- (fi|) are determined by successively solving the perturbed Einstein equa- 
tions with the same power of e ( |Hartle 1967| , |(Jhandrasekhar and Miller 1974] , |Quintana 1976| ). 

We consider the linear perturbations of the equilibrium fluid state. We use the Eulerian 
change of the pressure, density and 4-velocity, which are represented by Sp, Sp and 5U V . The 
perturbations of the gravitational field are neglected in the Cowling approximation. It is 
convenient in the following calculation to expand these functions in terms of the spherical 
harmonics YJ m , 



Sp 
Sp 

(p + p)5U r 
(p + p)SU e 



Irn 

^2$Plm(t,r)Yi m (9,(f; 

Irn 

e»/ 2 Y,Ri m (t,r)Y lm ( 



(5) 
(6) 
(7) 



W2 



Im 



E 

Irn 



V lm {t,r)d e Y lm (0A) ~ Ulm{ \ r) d^YU 

sin 6 
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(p + p)SU <p = e^ 2 J2[Vim(t,r)d^Y lm (e,<p) + U lm (t,r)smed e Y lm (e,<f ) )}, (9) 
i m 

su t = -nsuj,. (io) 

With these definitions, the pulsation equations are derived from the conservation laws of the 
perturbed energy-momentum, 

^„ = 0, (11) 

where 

5T; = (5p + 5p)£^" + [(p + p)<5t^tr + (p + p^tg + 5pS;. (12) 

Since we assume that the perturbation is adiabatic, the thermodynamic relation between 
the perturbed pressure and density can be written as 

5 P + Z-Vp=^-(5p + Z-Vp), (13) 
p + p 

where T is the adiabatic index and £ is the Lagrange displacement. 

We will solve eqs. (|TT|) and ( |i~3"D by the expansion of e. In a spherical star, the pertur- 
bations decouple into purely polar and purely axial modes for each I and m. A set of the 
functions Vi m = (Spim, $Pim, Rim, Vim) describes the polar mode, while the function Ai m = 
Ui m describes the axial mode. In the presence of the rotation, the mode will be mixed among 
the terms with different I, while the mode can still be specified by a single m (e.g., |Kojima 



1992| , [Kojima 1997|) . The coupled equations are schematically expressed as 

= [Aim] + £ X [Pl±lm] + £ 2 X [Ai m , Ai ±2 m) +■", (14) 
= [P lm ] +Sx [A l±lm ] +S 2 X [P lm , V l±2m ] + • • • , (15) 

where the symbol S means some functions of order e, and the square bracket formally 
represents the relation among the axial perturbation function A\ m , or the polar perturbation 
function Vi m appeared therein. We also assume that the time variation of the oscillation is 
slow and proportional to Q, i.e., d t ~ Q ~ 0(e). This is true in the r-mode oscillation, as 
will be confirmed soon. We look for the mode which is described by a single axial function 
with index (I, m) in the limit of e — > 0. That is, the polar part should vanish, while the axial 
part becomes finite. Hence, the perturbation functions are expanded as 

Am = + e 2 A^ + • ■ ■ , V lm = e(V$ + e 2 V% + ■■■). (16) 



Substituting these functions into eqs. (|14|) and ([Lq), and comparing the coefficients of e n (n = 
0, 1, 2), we have 

= [A$>], (17) 
= [eV^ + SxAi 1 ^ (18) 
= [e 2 A?2} + £ x [eVg\j + S 2 x [A^A^J = [e 2 A% + S 2 x A^\, (19) 
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where we have used Ai± 2m = and the relation fllS)) in the last part of (|I"9). Equation ( 1^) 
represents the axial oscillation at the lowest order, which can be specified by U^. Equation 
(|T9|) is the second-order form of it, and the term S 2 x can be regarded as the rotational 
corrections. We will show the explicit forms corresponding to eqs.(|T^)-([T9|) in subsequent 
sections. 



3. First-order solution 

The leading term of eq.flTTD is reduced to 

Uj$-im x u£>=0, (20) 

where 

x = WTT) w = WTT) {n - uj) > (21) 

and a dot denotes time derivative in the co-rotating frame, i.e., Ui m = (dt + imQ)Ui m . The 
evolution of the perturbation can be solved by the Laplace transformation, 

u(s,r)= I" Uj»(t,r)e-«dt. (22) 
Jo 

The Laplace transformation of eq.(|20|) is written as 

( s + im (Q - x )) u(s, r) - /« (r) = 0, (23) 

where describes the initial disturbance at t = 0. After solving u and using the inverse 
transformation, the solution in t-domain is easily constructed as 

U^(t,r) = J fi2(r) s + ir ^ n _ x) ds = f^(r)e-^-^H(t), (24) 

where H(t) is the Heaviside step function. We will consider t > region only, so that the 
function H(t) may well be omitted from now on. For the Newtonian star, m(Q~x) becomes 
a constant 



This is the r-mode frequency measured in the non-rotating frame ( Papaloizou and Pringlc 



1978 , Provost et al. (1981) , Saio 1982] ). In the relativistic stars, w is monotonically increasing 



function of r, vjq < w < wr. The possible range of the r-mode frequency is spread out. If one 
regards eq.(p4]) as the sum of the Fourier component e~ lat , then the spectrum is continuous 
in the range 

1 - — — --f mO < a < 1 - — — ma 26 

1(1 + 1) n J ~ ~ \ 1(1 + 1) n ) K ' 
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This result is the sameQ even if metric perturbations are considered ( Kojima 1998Q . The 
time-dependence is determined, whereas the radial dependence is arbitrary at this order. 
The function in eq.(^) is constrained by the equation of motions for the polar part, as 
will be shown in the subsequent sections. In this manner, the perturbation scheme (|r7|)-(|T§D 
is degenerate perturbation. 



4. Second-order equation 



The pressure and density perturbations are induced by the rotation, while the axial mode 
in the non-rotating stars is never coupled to them. We may assume that the perturbation 
can be specified by a single function l/j^ at the leading order, since general case can be 
described by the linear combination. The axial part induces the polar parts with index 
(/ ± l,m) according to the perturbation scheme fllSP . It is straightforward to solve Spi±i m , 
and 6pi±im from two components of eq.flTTD in terms of the first-order corrections and Uj^. 
The explicit results are 



where 



$Pl±lm 
8pi±lm 

s+ 
Q + 



2S+wU, 



(i) 

lm i 



-AS ± ( —^(e^wUl2)' + 2T z 



v 



I 



l + l 

IQ+, 



l + l 



e 
r v 



lm i 



I 



Q- 



m 



(l + l)Q- 

2 

Q- 



i 2 



\ (2Z-1)(2J + 1) 



(27) 
(28) 

(29) 
(30) 

(31) 



\ (2/ + l)(2/ + 3)' 

We here denote a derivative with respect to r by a prime. The lowest order form of eq. (|T3| 
is also decoupled into the equation of each l,m component, 



dp 



l±lm 



AC 2 e u e 



X Rl±lm 



-Q±u, 




with 



C 2 
A 



Po + Po ' 

Po _ Po_ 
Po + po Ypo 



(32) 

(33) 
(34) 



1 Note that this conclusion is derived within the lowest-order approximation of the perturbation scheme. 
The consistency in the higher order will affect it for barotropic case as will be shown in the following sections. 
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In eg. (|32|) , the displacement £ 2 represents the quadrupole deformation of the stationary star, 
and is related to the quantities of 0(e 2 ) as 



^2 



V 



2 2-v 



— [h 2 + -w r e 



(35) 



The region for A > is convectively unstable, while that for A < is stably stratified. We 
here assume A ^ in the following calculations, but will separately consider the case A = 
in Section 5. From eq . (|32|) , the function Ri±i m is solved by for the region A ^ 0. The 
function V\±\ m describing the horizontal motion is calculated from the remaining component 
of eq. (|TT|) . The explicit form is given by 



[l'{l' + l)V Vm ] v=l±1 



S ± (v 2 + 1(1 + lyne-U^) + T ± ( Vl + -r 2 we-uV) 
+ l(l + l)Q±v , 



(36) 



where 



v 2 



v 



4e -(3^+A)/2 



Av' 



2C 2 



v 



2e -(3v+\)/2 

3 6 r v(i)/ 



-_,(3u+X)/2 

Av' 



(r 2 zue 



I °lm 



r we 



-v/2\l 



u, 



2w 



lm 



2 

re 



ro + 2fi)--^--e^+ A )/ 2 
2 v ' 2 wr 2 



,(u+\)/2§2 
W 



U, 



(i) 



lm 



(37) 

Imi (38) 

(39) 



In eq.([J(J), we have used eq.(PUD to simplify it. 

These corrections in the polar functions affect the axial parts with indices (I ± 2, m) 
and (l,m). We are interested in the term with index (l,m) as the corrections to the leading 
equation. In deriving the axial equation with these corrections, the terms up to 0(e 3 ) in the 
background field also affect it. We include both corrections to the axial equation and have 
the following form, 

(40) 



= L/£> - im X U% + m 



(2) 



lm J 5 



where £ is the Sturm-Liouville differential operator defined by 



8c 3 we^ 2 -»\p +Po) 



r 2 e (X-u)/2 



?l 2 wU\ 



-Aczvj e 



Ac 2 w 2 e(- X -^' 2 



Po+Po 



,(A+3i/)/2 



(41) 



Av' 



w 



(wr 2 e u )'(po+Po) 
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G 



Az/r 2 



(42) 



^rV/ 2 Y 4c 2 2 4 



re 



3ci 

3m 2 
/(/ + !) 



-A/2 



- - — 6 - ^2 H ™2 H 

r A w r zu 



— + ^—6 + ^2 



5W 3 " 


_ m _ 




ti7 




w ) 



(43) 
(44) 



(I + l) n 

In order to solve eg. ([40"D , we introduce a complete set of functions as for the operator C, 

£k] + = 0. (45) 

where — k is the eigenvalue and the eigenfunction y K (r) is characterized by k, e.g., the number 
of the nodes. The eigenvalue is real number of 0(e 2 ), since L is the Hermitian operator of 
0(e 2 ). The eigenvalue problem is solved with appropriate boundary conditions. For example, 
the function should satisfy the regularity condition at the center and the Lagrangian pressure 
should vanish at the stellar surface. Certain initial date in eq.(p4|) can be decomposed 
by the set. We may restrict our consideration to a single function labeled by k, since the 
general case is described by discrete sum or integration over a certain range. By putting 
y K = imxfi^ lK} we have 

C[Ug>] = Ciimxf^e-^-^} = -imnxf^n^'^ , (46) 
where ( e -^-x)ty 

gives the first order correction and is neglected here. Using f^. K , we can 
integrate eq.(|40|) with respect to t, and have the function [7^(i, r) of 0(e 2 ) as 

Ug> = (imKxtfSb + £)e- im(Q - x)t , (47) 

where the function fj® of 0(e 2 ) is unknown at this order. The sum of the first and second 
order forms is approximated as 



U lm + U lm 



1 + imKxt)fi2,K + flm 



'f(l) , f( 2 ) 



-im(Q—x)t 



-im(Q,-(l+K)x)t 



(48) 
(49) 



(2) 

where we have exploited the freedom of / im to eliminate the unphysical growing term in 
eq. (|48"D . The value k originated from the fixing of the initial data becomes evident for large 
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t, since the accumulation of small effects from the higher order terms is no longer neglected. 
As a result, the frequency should be adjusted with the second-order correction to be a good 
approximation even for slightly large t, as in eq.(I9). This renormalization of the frequency 
is closely related to treating t as strained coordinate in the perturbation method. (See 
e.g. |Hinch 19910 In this way, the specification of the initial data at the leading order has 
influence on the second-order correction k. 



5. Oscillations in barotropic stars 

The structure of the neutron stars is almost approximated to be barotropic. The pul- 
sation equation is quite different from that of non-barotropic case, as in the Newtonian 
pulsation theory. The relation (j3^) for the case A = is replaced by 

OPl±lm = C dpi ±lm = —bpi± lm = — bpi m . (50) 

Po Vo 

In the last part of eq.(^), we have used hydrostatic equation of the non- rotating star. In 
this case, the function Ri±i m is never determined through eq.(|32|) unlike in A ^ case, but 
rather we have two restrictions to a single function . These conditions are never satisfied 
simultaneously unless for m = ±Z, in which one condition is trivial, 5p/_i m = 5p;_i m = due 
to Q_ = 0. The other condition for m = ±1 becomes 

«»r + (^-^<^)^^o. pi) 

Substituting the form (^4j) into this and neglecting the higher order term due to (e 



-im(n- x )ty 



we have the same differential equation for f^. The integration with respect to r results in 



W(+ r \ f(!) p -im{U- X )t 



Ulm r ) — flm:* e 



N (p + poy +l e- u (zue- u y l - 1)/2 \ e -*^"-W, (52) 

where A^o is a normalization constant. In this way, the function of the lowest order is 
determined. The corresponding the 3-velocity at t — > +0 is given by 

^ = 5v^ = N r l+ \we~^ 1 -^ 2 . (53) 



As shown by [Friedman and Morsink 1998| , the canonical energy of the perturbation is neg- 



ative for the Lagrange displacement and the solution therefore means unstable if the 
gravitational radiation reaction sets in. 



-10- 



We now specify Ri±i m of 0(e 2 ) to proceed to the pulsation equation with the second 
order corrections. The function X\±x m describing the radial motion is introduced as 

The function Xt±i m is arbitrary at this order, but we have two special cases. One is Xi±\ m = 
0, which is the limiting case for eq . fl32|) in a sense. The other is chosen so as to vanish 
the Lagrangian change of the pressure within the entire star. The condition corresponds to 
Xi±i m = 4:S±we x ^ u tj2^ I V '. For the first choice, the Lagrangian change of the pressure never 
vanishes at the surface unless po vanishes there. The second one is the usual way treated in 
the Newtonian r-mode, as shown by [Lindblom and lpser (1998)} With the second choice of 



Xi±im and eq.fl52"|), the pulsation equation for m = ±/ leads to 

- im X ui 2 2 - GUg* = 0, (55) 

where G is defined in eq.fl4"3"]). This equation is solved for llf® as in Section 4. The solution 
up to 0(e 2 ) is written as 



u lm ' Im — [J Im;* + J Im \ e 5 \ ou ) 

where the function of 0(e 2 ) is unknown at this order. This expression (|56|) is formally 
the same as in eq.([49|) for the non-barotropic case with k — G. The second-order correction in 

the frequency however depends on the position owing to the particular choice of the function 

f (i) 

J Im;* " 

We show the second-order rotational correction G for I = m = 2 mode in Fig.l. We 
adopt the polytropic stellar model with index n = 1. For the Newtonian star, G is a positive 
function, which monotonically increases from the center to the surface. The value ranges 
from G = 0.55(fi 2 R 3 /M) to G = 0.75(fi 2 i? 3 /M), where R and M are the radius and the 
mass for the non-rotating star. These values are rather smaller than that of incompressible 
case, G = 37 /27(Q 2 R 3 /M). (Only for the Newtonian incompressible case, the factor G 
is a constant. See Appendix.) The stellar deformation £ 2) which is the most important 
contribution to G, diminishes in the compressible fluid. As the star becomes relativistic, 
other relativistic factors become significant. As a result, the factor G is scaled down as 
a whole, and eventually becomes negative for some regions. In any cases, the frequencies 
satisfy the criterion of the radiation reaction instability, which implies < (1 + G)x/Q < 1- 
That is, retrograde in the rotating frame and prograde in the inertial frame. Therefore, the 
second-order correction never changes qualitative picture of the instability. 

In Fig. 2, we show the frequency range of the r-mode oscillations in the first-order ro- 
tational calculation. The upper and lower limits on the continuous spectrum in eq.(p6|) are 
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shown by two lines. The intermediate values between the two lines are allowed for a fixed 
model M/R. The frequency is a single value a^, given by eq.(|25|) in the Newtonian limit. 
The dragging effect relevant to the relativity broadens the allowed range as shown by eq.(E^). 



In Fig. 3, the allowed range is shown including the second-order corrections for the extreme 
case Q 2 = M/R 3 . Even in the Newtonian case, the oscillation frequency is not a single value, 
since the factor G depends on r, as seen in Fig.l. The allowed range of the frequency further 
increases with the relativistic factor. From this result, we expect that the r-mode frequency 
ranges from 0.8<7tv to 1.2<jn for a typical neutron star model with M/R ~ 0.2. 

We examine the effect on the spectrum of the gravitational waves emitted by the r-mode 
oscillations. We for simplicity neglect all relativistic corrections expect in the frequency, and 
estimate the spectrum by the Newtonian radiation theory. The dimensionless gravitational 
amplitude h(t) at infinity is determined by evaluating the time variation of the current 
multipole moment, h(t) ~ d l Si m /dt l flThorne 1980 ). The current multipole moment Si m for 



I = m = 2 mode and the velocity (|53|) is given by 

S 22 = N [ p (tue- u ) 1/2 e- 2l ^- x)t r 6 dr, (57) 



where the normalization N is also included the constant from the integration over angular 
parts. The Fourier component h(a) can be expressed as 

h(a) = f h(t)e +lat dt oc / p {we- u ) 1/2 (Q - xfK^ ~ 2(0 - x))r 6 dr. (58) 



The spectrum has a finite line breadth as shown in Fig.4. The spectrum is non-zero only for 
a = o"jv in the Newtonian treatment, but is broad in the relativistic one. The width at the 
half maximum is Act/un ~ 0.1. 



6. Discussion 

In this paper, we have calculated the r-mode oscillations in the relativistic rotating stars, 
neglecting the gravitational perturbations. The evolution can be described by oscillatory 
solutions which are neutral, i.e., never decay or grow in the absence of the dissipation. The 
oscillation is described not by a single frequency, but by the frequencies of a broad range, 
unlike in the Newtonian case. The reason is that the local rotation rate depends on the 
position due to the frame dragging effect even for uniform rotation. The r-mode oscillation 
frequency forms a continuous spectrum within a certain range. This property is distinguished 
from the well-known stellar oscillation modes such as the polar f-, p-modes, in which the 
spectrum of frequency is discrete. The r-mode frequencies lie in the unstable region for the 
gravitational radiation reaction instability, but it will be an important issue whether or not 
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the different spectrum of the unstable modes leads to different growth, e.g., in the non-linear 
regime. 
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A. Newtonian limit 



In this appendix, we will consider the r-mode oscillations in the Newtonian limit, in 
which 

(e u ,e\(p +p )/ Po ,zu/n) -1, (E/,m 2 ,fc 2 ,Wi,W 3 ) - 0. (Al) 



Equation (|40| ) reduces to the equation derived by Provost et al. ( 1981 )| , if the variables are 
matched. They solved the eigenvalue problem, assuming the form e -*( mfi - ff o(i+<Ti))t^ w here 
o"o of 0(e) is given by 2mfl/ (1(1 + 1)). The correction a\ of 0(e 2 ) was determined by solving 

in the case of A ^ 0, 



the eigenvalue problem for the operator C N = D N — (F N + Gn 

£n[v] + °\V = 0, 



with 



N 



G 



N 



Ac 3 n 2 Po 



-Ac^l 2 



Agpo 
' r 2 pp 



r 

,9 



+ 



Pa 

8c 2 n 2 r 
9 



Por 



+ 



Ag 
+ 2c x roL + 



'4 Cl tt 2 \ 

. A 9 ' ) 

2m 2 
1(1 + 1) 



a. 



where we have used the gravitational acceleration g = u'/2 and ellipticity a 
The second-order correction <j\ exactly corresponds to k in eq . (|45|) . 



(A2) 

(A3) 

(A4) 

(A5) 
-3&/(2r). 



As shown previously ( Provost et al. (1981)| , |Saio 198^ ), the eigen- value problem becomes 
singular^ for the barotropic case A = 0, since the second-order differential term vanishes. 
The second-order solution (1561) in the Newtonian limit reduces to 



U lm + u l- 



Im 



N oPo r l + fi 



Irn 



-i(mfl- (To(l+cri ))t 



(A6) 



2 Recently, 
stars. 



Lockitch and Friedman (1998) 



showed the resolution of the singularity in Newtonian isentropic 
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This equation is valid only for I = ±m, since the Newtonian counterpart of eq.([5(]) is never 
satisfied otherwise. When the correction is not a constant, the eigenvalue problem is 
ill-posed. The function G^ indeed depends on r for the compressible matter, so that |Provost| 
|et al. (1981) concluded no solution in this case. The exceptional case is the incompressible 
fluid, in which Gn is a constant since a = 5/4(Q 2 R 3 /M), g = Mr/R 3 . The correction in the 
frequency for I = m is 

Al n 2 R 3 21 

GN = -WTr^r + TTi a > (A7) 

which should be the same value calculated by [Provost et al. ( 1981 )| except a misprint in their 
expression. 
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Fig. 1. — The second-order correction in the frequency as a function of radius. The factor 
G is normalized by using the mass M, radius R and angular velocity Q. The solid line is for 
the Newtonian case M/R = 0.01. The dotted line represents the result for M/R = 0.1, and 
the dashed line for M/R = 0.25. 
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Fig. 2. — The frequency range with the relativistic factor M/R. Two lines denote the upper 
and lower limits of the frequency. The intermediate values between them are allowed for a 
fixed model M/R. The frequency is normalized by the Newtonian value <jn- 
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M/R 

Fig. 3. — The same as Fig.2, but the second-order correction with fl 2 = M/R 3 is included. 
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Fig. 4. — The Fourier transform of the gravitational wave amplitude of / = m = 2 mode for 
stellar model M/R = 0.2. The frequency is normalized by the Newtonian value <jn, and the 
amplitude h(a) is normalized by the maximum. 



